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Abstract 

Let G = {VG, AG) be a directed graph with a set S C VG of terminals and 
nonnegative integer arc capacities c. A feasible multiflow is a nonnegative real 
function F(P) of "flows" on paths P connecting distinct terminals such that the 
sum of flows through each arc a does not exceed c(a). Given /i: S x S ^ M+, 
the fi-value of F is Y2p F{P)fi{sp,tp), where sp and tp are the start and end 
vertices of a path P, respectively. 

Using a sophisticated topological approach, Hirai and Koichi showed that 
the maximum //-value multiflow problem has an integer optimal solution when 
fi is the distance generated by subtrees of a weighted directed tree and (G, S, c) 
satisfies certain Eulerian conditions. 

We give a combinatorial proof of that result and devise a strongly polynomial 
combinatorial algorithm. 

Keywords: directed multiflow, tree-induced distance, strongly polynomial algo- 
rithm 

1 Introduction 

1.1 Multiflows in directed networks 

We use standard terminology of graph and flow theory. For a digraph G, the sets 
of its vertices and arcs are denoted by VG and AG, respectively. A similar notation 
is used for paths, cycles, and etc. For X G VG, the set of arcs of G entering (resp. 
leaving) X is denoted by Sq{X) (resp. 6^^{X)). When X = {v}, we write Sq{v) for 
Sq{{v}). When G is clear from the context, it is omitted from notation. Also for a 
set A and a singleton a, we will write A — a for A \ {a}, and AU a for AU {a}. 
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A directed network is a triple M = {G, S, c) consisting of a digraph G, a set 
S C VG of terminals, and integer arc capacities c: AG — )■ Vertices in l^G — S 
are called inner. A directed path in G is called an S-path if its endvertices are distinct 
elements of S. A multiflow F is a function assigning a nonnegative real number, or 
flow, to each S-path. A multiflow F is called feasible if for each arc a G AG, the sum 
of flows assigned to S-paths going through a does not exceed c(a). The (total) value 
of F is the sum of flows over all S-paths P: 

(1.1) val(F):=^^F(P). 

Sometimes (e.g., in |IKN98] ) such multiflows are called free to emphasize the fact 
that any pair of terminals is allowed to be connected by nonzero flows. 
The following maximum (fractional) multiflow problem is well known: 

(MF) Given a directed network M = {G, S, c), find a feasible multiflow F of maxi- 
mum value. 

The problem in which one is asked for maximizing among the integer multiflows is 
denoted by IMF. 

For general directed networks J\f, problem IMF is NP-hard already for 15*1 = 2 
|FHW80] . Tractable cases have been revealed for networks obeying a certain con- 
servation property. More precisely, c (or A/") is called Eulerian at a vertex v if 
c((5™(f)) = c(5°"*(f)). (For a function / : A — )■ M and a subset A' C A, we write 
f{A') for ^(/(a): a G A').) When c is Eulerian at all inner vertices (resp. at all 
vertices), the network A^ is called inner (resp. totally) Eulerian. By a cut in G we 
mean a pair of nonempty subsets {X,X), where X C VG and X := VG — X. It is 
called an (5*1, S2)-cut if Si C X and 5*2 C X. When it is not confusing, we may refer 
to the arc sets 5'q^{X) and 5q{X) as cuts as well. 

The following result signifies the importance of inner Eulerian networks: 

Theorem 1.1 (Lomonosov (unpublished, 1978), Frank |Fr89j ) Let M = 

{G, S, c) be an inner Eulerian directed network. Then there exists an integer max- 
imum feasible multiflow F in M . It satisfies 

val(F) = 5^^^^c(5-*(X,)), 

where for each t E S , {Xf, Xf) is a minimum capacity {t, S — t)-cut in M . 

Therefore, MF and IMF have the same optimal value for an inner Eulerian 
network, and this value can be found in strongly polynomial time (by computing a 
minimum (t, S — t)-cut for each t E S). Ibaraki, Karzanov, and Nagamochi |IKN98| 
devised a "divide-and-conquer" method that computes an integer maximum multiflow 
in such anetwork in 0((MF(n, m)+mn)-log ISI+mn^) time. (Hereinafter n := \VG\, 
m := \AG\, and MF{n' ,m') denotes the complexity of a max-flow computation in a 
directed network with n' vertices and m' arcs.) The latter complexity was improved 
to 0{{MF{n,m) + mn\og{n^/m)) -logl^l) in |BK07j . 
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1.2 Weighted multiflows 

A generalization of MF involves weights between terminals. More precisely, given a 
weighting fi: S x S R+, the fi-value of a multiflow F is 

(1.2) val(F,/i) := /i(sp, tp)F(P), 

where the sum is over all S-paths P, and sp and tp denote the start and end vertices 
of P, respectively. We may assume that fi{s, s) = for all s G S*. 

Replacing (II. ip by fll.2p . we obtain the weighted counterpart of MF: 

(/x-MF) Given M and fi as above, find a feasible multiflow F of maximum fi-value. 

The integer strengthening of yU-MF is denoted by /i-IMF. When fi{s, t) = 1 for all 
s^t, n-MF turns into MF, and /x-IMF into IMF. 

1.3 Tree-induced weights 

It has been shown that problem yU-IMF has a rather wide spectrum of tractable cases. 
The simplest case is 5* = {s, t}, n{s, t) = 1 and n{t, s) = 0; then /i-IMF becomes the 
standard maximum flow problem with arbitrary integer capacities. A representative 
well-solvable class has been found in connection with the so-called directed multiflow 
locking problem: 

(DMLP) Given a directed network M = [G, S, c) and a collection C C 2'^, find a 
feasible multiflow F in M that locks simultaneously all members of C. 

Here F is said to lock a subset A (Z S ii the sum of values F{P) over the S-paths P 
going from A to 5* — A is maximum possible, i.e. it is equal to the minimum capacity 
of an (A, S — y4)-cut in M . A collection C C 2'^ is called lockable if DMLP has a 
solution for all (G, c) (with S fixed). Important facts are given in the following 

Theorem 1.2 ((iKNp]) C CI 2^ is lockable if and only if C is cross- free, i.e. 
for any A,B E C, at least one of the following holds: ACB,BCA, Ar\B = (l}, 
AUB = S . Moreover, ifC is cross-free and is inner Eulerian, then DMLP has an 
integer solution. Such a solution can be found in 0{{MF{n, m) + mn) ■ log 15*1 +mn^) 
time. 

This gives rise to the following tractable cases of yU-IMF. Given C C 2"^, take an 
arbitrary function i : C ^ ]R_|_. For s,t E S, define 

(1.3) /i^(s, t) := ■■ AeC, seA^t). 

Suppose that C is cross-free and Af is inner Eulerian, and let F be an integer solution 
to DMLP (existing by Theorem II. 2p . Then F is simultaneously an optimal solution 
to /i^-IMF for every £ : C — > R+; this can be easily concluded from the fact that F 
saturates minimum capacity {A, S — y4)-cuts in J\f for all A E C. 

A cross-free collection C can be represented by use of a directed tree T = {VT, AT) 
(a digraph whose underlying undirected graph is a tree); namely: 
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(1.4) there is a bijection /3 : C — AT and a map 7 : 5 — )■ VT such that for each arc 
a = {u,v) G AT and for A := /3~^(a), the set of terminals s G 5* whose image 
7(3) occurs in the component of T — a containing u is exactly A. 

Under this correspondence, we may interpret £ : C — )■ as a length function on the 
arcs of T, keeping the same notation: i{a) := i{A) for a G AT and A := /3~^(a). 
These arc lengths induce distances d = di on VT in a natural way: 

(1.5) for x,y & VT, define d{x, y) to be the sum of ^-lengths of forward arcs in the 
simple path from x io y mT . 

(This path may contain both forward and backward arcs. If there are no forward 
arcs, we have d{x,y) = 0.) One can see that fi = fi£ figured in fll.3p satisfy 

(1.6) fi{s,t) = d,{-f{s),j{t)) s,teS. 

Note that such a fi satisfies the triangle inequalities fi{s,t) + fi{t,u) > fi{s,u) for 
all s,t,u G S, i.e. is a directed metric (space). In light of f ll.5p -( ITT6|) . fj, is called a 
tree-induced directed metric. 

Generalizing the above-mentioned integrality results, Hirai and Koichi |HK10| 
considered arbitrary weight (or distance) functions fi : S x S M+ and gave an ex- 
haustive analysis of the integrality and "unbounded fractionality" behavior of problem 
/i-MF in terms of yU, for both integer and Eulerian cases. 

More precisely, let T be a directed tree with nonnegative arc lengths i. Given 
a network JV as before, suppose that each terminal s G S* is associated with some 
subtree (a weakly connected subgraph) of T. For s, t G 5*, define fi{s, t) to be the 
distance from Tg to Tt, i.e. 

/i(s, t) := mm{d{u, v) : u G Tg, v ^ Tt}. 

(Note that such a fi need not satisfy triangle inequalities. When each is a single 
vertex, /i is specified as in (11. 6p .) 

A weight function fi on S x S that can be obtained in this way is called a tree- 
induced (directed) distance, and an appropriate tuple TZ = (T, i, {Tg}) is called a (tree) 
realization of fi. 

For such an TZ, we distinguish between three sorts of terminals. We call s G S* 
simple if Tg consists of a single vertex, linear if Ts is a directed path, and complex 
otherwise. If the whole tree T forms a directed path, then TZ itself is called linear. 

Hirai and Koichi obtained the following results. 

Theorem 1.3 ( |HK10| ) LetM = {G, S, c) be a directed network and let 11: SxS — )■ 

M_i_ be a tree-induced distance with a realization TZ. 

(i) IfTZ is linear, then fi-MF has an integer optimal solution. 
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(ii) If c is Eulerian at all inner vertices and all complex terminals, then fi-MF has 
an integer optimal solution. 



Note that the proof of (i) in Theorem 11.31 given in |HK10] is constructive; it 
reduces problem /x-MF in this case to finding a certain min-cost circulation. The 
proof in [HKIOJ is much more involved; it employs a topological approach based on a 
concept of tight spans of directed distance spaces introduced in that paper. (Another 
nice result in |HK10) relying on the directed tight span approach asserts that if a 
distance /i is not tree-induced, then /i-MF has unbounded fractionality in the totally 
Eulerian case, i.e. there is no positive integer k such that /i-MF admits a ^-integer 
solution for every totally Eulerian network and this yU.) 

In this paper we devise an efficient combinatorial algorithm that constructs an 
integer optimal solution to problem /i-MF under the conditions as in (ii) of The- 
orem II. 3t this yields an alternative (and relatively simple) proof of assertion (ii). 
Our method extends the divide-and-conquer approach of |IKN98) : it is described in 
Section [21 The algorithm runs in 0{{MF{n,m) + mn\og{n? / m)) ■ log \S\) time. 

2 Algorithm 

Let M = {G, S, c) and /i obey the conditions in Theorem II. 3( ii). 

The following convention will allow us to slightly simplify the description of our 
algorithm (without loss of generality). In a tree realization TZ = {T,i, {Ts}) of a 
distance function /i, let us think of T as an undirected tree with edge set ET, and 
assume that each edge e = uv generates two oppositely directed arcs: one going from 
u to V, and the other from f to m (yielding a "directed quasi-tree") . For a = {u, v), the 
opposite arc [v, u) is denoted by a. The length function £ is given on the corresponding 
arc set, denoted by AT as before, and for x,y E VT, the distance d{x,y) is defined 
to be the ^-length of the corresponding directed path from x to y. Accordingly, a 
terminal s is linear if the subtree Tg is a simple undirected path and one of the two 
directed paths behind Tg has zero ^-length. 

Sometimes, to ensure the desired efficiency of the method, we will be forced to 
treat some linear terminals as complex ones (which will never be confusing). 

2.1 Initial reductions 

Let TZ = (T, i, {Ts}) be a tree realization of /i. A pre-processing stage of the algorithm 
applies certain reductions to TZ (called initial reductions) . 

Choose a linear terminal s G S* (if exists), i.e. Tg is a path in T connecting some 
vertices ti and t2, and one of the two directed paths behind T^, from ^2 to ti say, has 
zero ^-length. This implies that 

(2.1) fj,{x, s) = d{Tx,ti) and fj.{s,x) = d(t2,Tx) for all a; G 5. 
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We replace s by a pair of simple terminals as follows. Add to G new terminals 
Si and S2 and arcs (s, Si) and {s2,s). The capacities of these arcs are chosen to be 
sufficiently large and to make the network Eulerian at s. Denote the resulting digraph 
by G' and the network by A/"' := (G", S', c'), where S' := {S — s)U {si, 32}- We modify 
TZ into the tuple TZ' with the same tree T by setting T^^ := {ti} and := {^2}- This 
gives new distance /x': S" x S" — )■ M+. 

We claim that the two problems: /i-MF with (A/',7^) and /i'-MF with (A/'^7^'), 
are essentially equivalent. Indeed, /i'(s2,Si) = (since the ^-length of the directed 
path from t2 to ti in T is zero). Therefore, one may consider only those multiflows 
in A/"' that are zero on all S2-S1 paths. Any other 5"-path P' in A/"' has a natural 
image (an S-path) P in J\f. Namely, if P' neither starts at S2 nor ends at si, then 
P = P'. If P' starts at S2 (resp. ends at Si), then P is its maximal subpath from 
s (resp. to s). This gives a one-to-one correspondence between the S"-paths P' in 
A/"' (excluding S2-S1 ones) and the S'-paths P in A^, and by (12.11) . the transformation 
preserves distances: fi'{sp',tp') = fi{sp,tp). We reset J\f := Af' and 71 := 71'. 

Making a sequence of similar reductions, we obtain a situation when 

(CI) Any terminal in A/" is either simple or complex. 

Four more sorts of reductions are applied to ensure the following additional properties: 

(C2) Each leaf in T corresponds to some (possibly multiple) simple terminal. (For 
otherwise the leaf can be removed from T.) 

(C3) No inner vertex f of T corresponds to a simple terminal. (For otherwise one 
can add to T a new vertex v' and edge vv' with zero ^-length of both arcs 
{v,v'), {v',v), and replace the subtree {v} by {v'} in the realization of /i.) 

(C4) Each inner vertex of T has degree at most 3. (This can be achieved by splitting 
inner vertices of bigger degrees in T and by adding additional edges with zero 
£- length of arcs in both directions.) 

(C5) T has 0(1 vertices. 

To provide (C5), note that T has OdS*!) leaves (by (C2)), and hence it has OdS*!) 
vertices of degree 3. Consider a vertex v of degree 2 in T, and let e = uv and e' = vw 
be its incident edges. If for any s E S, the subtree contains either none or both 
of e, e', then we can merge e, e' into one edge uw (adding up the corresponding arc 
lengths), obtaining a realization of fi with a smaller tree size. Otherwise f is a leaf 
of some Tg. Obviously, the number of leaves of Tg does not exceed that of T, so it is 
estimated as OdS*!). This gives (C5). 

2.2 Optimality certificate 

Here we establish a sufficient condition that implies optimality of a given multifiow. 
We need some additional terminology and notation. A feasible multifiow F in A/" is 
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said to saturate a cut {X,X) in G if each S-path P with F{P) > meets U 
at most once, and 

V F(P) = c(e) for each arc e G 

Definition. For an arc a = {u, v) of T, define Ua to be the set of pairs [s, t) in S 
such that fi{s,t) "feels" i{a), i.e. a belongs to a minimal directed path that starts in 
Ts and ends in Tt. (Then n{s,t) is just the ^-length of such path.) 

For a multifiow F and a set 11 C 5* x 5, let -F[n] be the "restriction" of F relative 
to n. More precisely, for an S'-path P in G, we define 



F[n](P) : = 



F{P) ii{sp,tp)eu, 

otherwise. 



A cut (X, X) in G is called U-separating if s G X ^ t holds for each (s, t) G 11. 
Lemma 2.1 Let F be a feasible multifiow in M . Suppose that 

(2.2) there exists a collection {(Xa,Xa): a G AT} of cuts in G such that for each 
a G AT, (Xa,Xa) is a Ua- separating cut saturated by F. 

Then F is an optimal solution to fi-MF. 

Proof For s,t E S, let f{s,t) denote the sum of flows (by F) over the paths from s 
to t in G. Then 



Also 



/i(s,t) = y i{a). 



It follows that 

(2.3) --KF,,) = j:^^^^e{a)[j:^^^^^^^J{s,t) 

Consider an arc a G AT. Since {Xa,Xa) is a Ila-separating cut, we have 

(2.4) E„,„„./(».*)Sc(«-(^J). 
Then (Q and (EH) give 

(2.5) val(F,/i)<5^ £(a)c(5-HX.)). 

Since each cut {X^, Xa) is saturated by F, inequality 02. 4p turns into equality, and 
so does (12. 5p . Thus, val(F, /i) is maximum, and the lemma follows. □ 

Given a problem instance (A/", TZ), our algorithm will construct an integer multiflow 
F that possesses property (12. 2p . and therefore F is optimal by Lemma 12.11 Note 
that (12. 2p does not involve the lengths i of arcs in T, so F is optimal simultaneously 
for all distances /i induced by arbitrary i (when T and {Tg} are fixed). 
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2.3 Partitioning step 

The core of the algorithm consists in the following recursive procedure that divides 
the current instance (A^, TZ) into a pair of smaller ones. 

Suppose that T contains an edge e = V1V2 such that neither Vi nor V2 is a leaf. 
Let a := (fi,f2)- Deletion of e splits T into subtrees Ti and T2 with Vi G VTi and 
V2 G VT2. Define Si (resp. S2) to be the set of terminals s & S such that is 
entirely contained in VTi (resp. in VT2). Then fl 6*2 = and each terminal in 
S — {Si U 5*2) is complex (by properties (C1),(C3)). Hence JV is Eulerian at each 
vertex in VG — [Si U 6*2). Also from the definition of Ua it follows that 

Ua = SixS2 and Ua = S2xSi. 

Compute an (5*1, S'2)-cut {Xi,X2) of minimum capacity c(5°^*(Xi)) in G. Then 
{Xi,X2) is Ha-separating and {X2,Xi) is Ha-separating. The Eulerianess implies 

c{6^-\Xi)) - c(5-(X0) = J2,, - 

Hence the capacity c{6°^^{X2)) = c(5™(Xi)) is minimum among all (5*2, S'i)-cuts in G 
as well. 

We construct two new instances (A/i,/ii) and 1/^2, 1^2) in a natural way. More 
precisely, set A/i := (Gi, S"^, Ci), where Gi is obtained from G by contracting X2 into 
a new vertex Z2 (and deleting the loops if appeared), Ci is the restriction of c to the 
arc set of Gi, and S'l := (S* fl Xi) U {-22}- The distance fii is induced by the tree 
obtained from T by contracting the subtree T2 into t'2; the arc lengths in T[ are same 
as in T. (In fact, these lengths are ignored by the algorithm and they are needed 
only for our analysis.) Terminals s G 5* fl Xi are now realized by the subtrees of T[ 
obtained by restricting the subtrees Ts in TZ to T[. The terminal Z2 is realized by 
{^2}. Let TZi denote the resulting realization of fii. 

The construction of Af2 = {G2, 5*2, C2), fi2, is symmetric (by swapping 1^2). 

The algorithm recursively constructs integer optimal multifiows Fi and F2 for 
(AfijTZi) and (A/'2,7^2), respectively. The following property easily follows from the 
minimality of (Xi,X2) and (X2,Xi): 

(2.6) the multifiow Fi saturates the cuts 6q {Z2) and 6q^^{z2); similarly, F2 saturates 
<5j,-(^i)and5-;(^i)- 

This property enables us to "glue" (or "aggregate") Fi and F2 into an integer mul- 
tifiow F in A/" which saturates both cuts (Xi,X2) and (X2,Xi) (being Ha-separating 
and Ha -separating cuts, respectively). These cuts together with the preimages in G 
of corresponding saturated cuts for Fi and F2 give a collection of saturated cuts for 
F as required in (12.21) . yielding the optimality of F by Lemma [2.11 

The above partitioning step reduces the current problem instance to a pair of 
smaller ones (in particular, the tree sizes strictly decrease). One easily checks that 
conditions (C1)-(C5) (see Section I2TT]) are maintained. Note that for i G {1,2}, if s 
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is a complex terminal in S such that the image of Tg in T/ is different from {f3_j}, 
then we should keep regarding s as a complex terminal in A/i (even if this image 
is a nontrivial (undirected) path having zero ^-length in one direction). This is not 
confusing since the network continues to be Eulerian at s. 

The recursion process with a current T stops when each edge in it is incident to 
a leaf. Since each inner vertex of T has degree 3 (by (C4)), only two cases of T are 
possible: 

(i) VT consists of two vertices Vi and V2', 

(ii) VT consists of one inner vertex vq and three leaves fi,f2,f3. 

Case (i) is considered in Subsection 12. 4^ and case (ii) in Subsection 12.51 

2.4 Basic step: two vertices 

Let e = V1V2 be the only edge of T. Note that the vertices f 1 and f 2 may correspond 
to many terminals in S. Let terminals Si, . . . ,Sp (resp. ti, . . . ,tq) be realized in TZ by 
{vi} (resp. {^2})- Also there may exist a terminals s realized by the whole tree T; 
but such an s may be ignored since fi{s, t) = fi{t, s) = for any t G 5*. 

Let S' := {si, . . . , Sp} and T' := {ti, . . . , tg}. Construct an integer maximum S'- 
T' flow, i.e. a function /: AG — ?■ Z_|_ with val(/) := ^ (divj(s) : s G S") maximum 
subject to /(a) < c(a) for each a G AG and divj(f ) = for each v G VG — {S' U T'). 
Here divj(f) denotes the divergence f{S°^^{v)) — f{6™{v)). Then / saturates some 
(S",T')-cut (X, X). Since capacities c are Eulerian at all inner vertices, g := c — f is 
a T'-S' flow. This implies that val{g) = c{X, X) and that g saturates the reversed 
cut (X,X). 

We construct F by combining path decompositions of / and g. Let a := (fi,f2)- 
Then Ua = S' x T', and {X,X) is a Ila-separating cut. The multifiow -Fpa] corre- 
sponds to /, and therefore it saturates {X,X). Similarly, IIq: = T' x 5", {X,X) is a 
Ila -separating cut, the multifiow -Fpa] corresponds to g, and therefore it saturates 
(X,X). This gives (O for F. 

2.5 Basic step: three leaves 

This case is less trivial. Here FT consists of three edges = ViVo, i = 1,2,3. We 
denote the arc (fi,fo) by Oj. 

Let us call terminals s, s' in the current network = {G, S, c) similar if they are 
realized by the same subtree of T; clearly fJ,{s,p) = fj,{s',p) and /i(p, s) = /i(p, s') for 
any p E S. Suppose that there are similar simple terminals s, s'. They correspond to 
the singleton {vi} for some i G {1,2,3} (in view of (C3)). The fact that Vi is a leaf 
of T provides the triangle inequality fi{p, s) + fi{s, q) > fi{p, q) for any p,q E S, and 
similarly for s'. Due to this, we can identify s, s' in G into one terminal (corresponding 
to {vi}) without affecting the problem in essence. 
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Thus, we may assume that for each i = 1,2,3, there is exactly one terminal. 
Si say, corresponding to {fj}. Let S' := {si, 52,53}. Note that each terminal s G 
5* — S" is (regarded as) complex, and c is Eulerian at s. We partition S — S' into 
subsets 5*1, 5*2, 5*3, Su, S13, S23, where Si (resp. Sij) consists of the (similar) terminals 
corresponding to the subtree of T induced by the edge Cj (resp. by the pair {cj, ej}). 

Suppose we ignore the terminals in S" — S', by considering the network A/"' := 
(G, S', c). This network is inner Eulerian since Af is Eulerian within S — S'. Using 
the algorithm from |BK07] . we find an optimal multiflow F to problem IMF for A/"' 
with unit distance for each pair (sj, Sj), i 7^ j. Also for z = 1, 2, 3, we find a minimum 
capacity (sj, S'—Si)-cut (Xj, Xj) in A^'. They can be chosen so that the sets Xi, X2, X3 
are pairwise disjoint. Also one may assume that each path P with F{P) > is simple 
and has no intermediate vertex in 5". Then F yields a solution to (A/"', /i'), where fi' is 
the restriction of utoS'x S'. Since A/"' is inner Eulerian and in view of Theorem 1 1.1^ 
F saturates both cuts (Xj,Xj) and (Xj,Xj) for each i. Associating such cuts to the 
arcs ai,ai results in (12. 2p . Then F is optimal by Lemma [2.11 

Next we return to A/" as before. The above multiflow F need not be optimal 
for {M 1 ^) since cuts (Xj,Xj) may not be na--separating for some i. Our aim is to 
improve F, Xi , X2 , X3 so as to ensure (12. 2 p . 

More precisely, we are looking for subsets X[ C Xj, i = 1, 2, 3, and a multiflow F' 
such that: 

(2.7) (i) SiCVG - {Xj U Xfc) and Sij CVG-Xk for any distinct k; 

(ii) for i = 1,2,3, the cuts (X'-,X'-) and {X'-,X^) are saturated by F'; 

(iii) each path P with F'{P) > connects either Si and Sj, or Si and Sj, or Sj 
and Sjk, where k are distinct. 

By ([22D(i), the cut (X;,X^) is n„^-separating, i = 1,2,3. In their turn (122D(ii),(iii) 
imply that F'paJ saturates {X'-,X-). Then F' is optimal by Lemma [2.11 

We construct the desired X^' and F' as follows. For i = 1,2, 3, let Qi denote the set 
of terminals s that violate fl2.7p (i) w.r.t. X,, i.e. s G Xj but s ^ {sj} U 5*^ U Sij U Sik 
(where {i,j,k} = {1,2,3}). If = then (Xj,Xj) is already na.-separating, in 
which case we set X^' := Xj. 

Let Qi 7^ 0. We construct the digraph Gi from G by contracting VG — Xi into 
a new terminal Zi. Arc capacities in Gi are induced by those in G (and are denoted 
by c as before). This gives the network A/j := {Gi, {si,Zi} U Qi, c) which is Eulerian 
at all vertices except, possibly, for Si and Zi. The current multiflow F in G induces 
a multiflow Fi in Gi consisting of weighted Si-Zi and Zj-Sj paths. Since F saturates 
S'SiXi) and 6^\Xi), the multiflow Fi saturates 6'S^{zi) and 6^f{zi). 

Now we find in Gi a maximum integer flow Qi from the source Sj to the set of sinks 
Qi U Zi. Moreover, among such flows we choose one maximizing —dwg.{zi). (This is 
done by standard flow techniques: take a maximum Si-Zi flow (e.g. by extracting the 
subflow in Fi formed by Sj-Zj paths), then switch to the residual network and augment 
the current flow to get a maximum Si-{Qi U Zi) flow.) The flow Qi is decomposed into 
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a collection of weighted Si-Zi paths, denoted by Zi), a collection of weighted Sj-t 
paths for t G Qi, denoted by ^i{si,t), ignoring possible cycles. By the construction, 
gi saturates the trivial cut ^^'^ some (sj, Qi U 2;j)-cut 

It remains to construct flows on paths going in the opposite direction, i.e. entering 
Si. Define the function hi := c — Qi on AGi. It is Eulerian at all vertices in VGi — 
{{si, Zi} U Qi). Also divhiizi) > (since Qi saturates 6Q.{^i)) div/i.(t) > for 
all t & Qi (since c is Eulerian at t). We decompose hi into a collection of weighted 
ZiSi paths, denoted by hi{zi,Si)), and a collection of weighted t-Si paths for t G Qi, 
denoted by hi(t,Si), ignoring possible cycles. These paths saturate 6q.{XI) (since 
hiiS^J^iX^)) = 0) and 5S|*(^,Hsmce h,{S'SSz)) = 0). 

The collections ^i{-) and hi{-) constitute a multifiow F/ that replaces the "restric- 
tion" Fi of F on Gi. Making such "replacements" for i = 1,2,3 (and using the fact 
that Xi, X2, are disjoint), we obtain an integer multifiow F' which along with 
X[,X2,X^ as above satisfies ( 12. 7p . Hence F' is optimal. 

2.6 Complexity 

In this final section we describe an efficient implementation of our algorithm and 
estimate its complexity. Current multifiows in the process are stored as collections 
of point-to-point flows. Namely, an integer multifiow F in a network with terminals 
S is maintained as a collection \ s,t & S,s ^ t}, where f^t is an integer s-t flow 
(called an s-t component of F). 

Let ip{n, m, k) denote the complexity of the algorithm applied to an instance with 
n vertices and m arcs of G, and k leaves of T. 

The case k = 2 was studied in Subsection 12.41 The algorithm involves a single 
max-flow computation and two flow decompositions. Hence 

(2.8) Lp{n,m,2) = 0{MF{n,m) +mn), 

where MF{n', m') denotes the complexity of a max-flow algorithm in a network with 
n' vertices and m' arcs. 

The case k = 3 was considered in Subsection 12.51 It reduces to solving a three- 
terminal version of the unweighted directed IMF problem followed by 0(1) max-flow 
computations and decompositions. With the help of the algorithm from f BK07] the 
three-terminal multifiow problem is solved in 0{MF{n,m) + mnlog(n^/m)) time. 
Therefore, 

(2.9) (p(n,m,3) = 0{MF{n,m) + mn\og{n^ /m)). 

For /c > 4, we apply the partitioning operation from Subsection 12.31 Computing 
a minimum cut dividing the current instance M into Mi and N'2 takes 0{MF{n, m)) 
time. The aggregation takes the S-Z2 components of Fi (for s G SClXi) and the ZiS 
components of F2 (for s G fl X2), combines them into an (5* fl Xi)-{S fl X2) flow 
and decomposes it into a collection of flows for all source-sink pairs. The algorithm 
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similarly handles the Z2-s components of Fx (for s G STlXi) and the s-Z\ components 
of F2 (for s e 5" n X2). Finally it adds remaining components of Fx and F2, thus 
forming an integer optimal multiflow in M . In total the aggregation operations take 
Oimn) time, hence 

(2.10) V9(n, m, A;) = ^{n\, mi, fci) + ^{n^-, "^2, ^2) + 0{MF{n, m) + mn), 

where ( the size parameters for Mi. 

Since degrees of inner nodes of T are 3 by (C3), there exists (and can be found 
in 0{k) time) a partitioning edge in T that yields fci,/i;2 ^ 2/c/3 + 1. Thus, 
the height of the recursion tree is at most 0(\ogS). Also rii + n2 = n + 2 and 
Tni,m2 < m. Assuming that MF{n,m) obeys some technical conditions (e.g., sat- 
isfying MF{n,m) = 0(mn log (n^/m)) as in the algorithm of Goldberg and Tarjan 
|(;T88] ). one can show by induction that (^M), (O, and fl2:T0|) imply 

ip{n, m, k) = 0{{MF{n, m) + mn) log k + mn log(n^/m)) 

(applying reasonings similar to those in |IKN98) ). By spending additional 
0{mn log 15*1) time, one can convert the resulting integer optimal multiflow into path- 
packing form, as exlained in {BK07] . In total, the algorithm takes 0{{MF{n,m) + 
mn) log 15*1 + mn log(n^/m)) time, as declared. 
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